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FOREWORD 


This report is prepared by the Lockheeed-Georgia Company under Con- 
tract NASI-159^9. "Advanced Composite Structural Design Technology for 
Commercial Transport Aircraft," and describes the theory and capability of 
a computer program prepared for the analysis and sizing of stiffened com- 
posites panels. This work was performed under Task Assignment No. 5 of the 
contract. The program is sponsored by the National Aeronautics and Space 
Administration, Langley Research Center (NASA/LaRC). Dr. James H. Starnes 
is the Project Engineer for NASA/LaRC. John N. Dickson is the Program 
Manager for the Lockheed-Georgia Company. 

In addition to the authors the following Lockheed specialist/consult- 
ants made major contributions to the material presented. 

L.W. Liu Programming 

Dr. J.T.S. Wang (Georgia Tech.) Analysis 
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POSTOP: Postbuckled Open -ST 1 f fener Optimum Panels - Theory and Capability 

John N. Dickson and Sherrill B. Biggers 
Lockheed-Georgia Company 
Marietta, Georgia 


SUMMARY 

A computer program, POSTOP, for the analysis or sizing or stiffened 
panels is described. Buckling resistant or postbuckled panels subject to 
loading conditions typical of aerospace structures may be treated. Com- 
posite materials may be used. Analytical routines are included to compute 
panel stiffnesses, strains, local and panel buckling loads, and skin/ 
stiffener interface stresses. Multiple load cases may be defined. Sizing 
i3 performed by the COPES/CONMIN optimization routines. This report 
describes the analytical and sizing procedures used in the program. 

INTRODUCTION 


Stiffened panels are widely used in aerospace structures. Applica- 
tions include fuselages, lifting and control surfaces, spar webs, bulkheads 
and floors among others. The skin or web is typically required to remain 
unbuckled under some loading conditions but allowed to enter the postbuckl- 
ing regime under other more severe loading conditions. For example, metal 
transport fuselages are normally allowed to buckle at ultimate loads but 
are required to be unbuckled at loads corresponding to 1g. level flight. 
The preliminary design of minimum-weight stiffened panels subject to multi- 
ple loading conditions and various requirements associated with each load- 
ing condition can become a formidable problem. If composite materials are 
to be considered, the number of design variables and the size of the opti- 
mization problem, as well as the complexity of the required analyses, may 
increase substantially. 
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As an aid to the preliminary design of stiffened composite panels that 
may be loaded in the postbucklirtg regime, an analysis and Sizing code, 
POSTOP, has been developed. A comprehensive set of analysis routines has 
been coupled to the COPES/CONHIN (Reference 1) optimization program to pro- 
duce this sizing code. POSTOP is intended for the preliminary design of 
metal or composite panels with open -sect ion stiffeners, subjected to com- 
bined biaxial compression (or tension) , shear and normal pressure loading. 
Longitudinal compression, however, is assumed to be the dominant loading. 
Temperature, Initial bow eccentricity and load eccentricity effects are 
included . The panel geometry is assumed to be repetitive over Several bays 
in the longitudinal (stiffener) direction as well as in the transverse 
direction. These restrictions were imposed to allow the development Of* 
analysis routines with the efficiency required in a practical Sizing eod§.- 
The resulting program is applicable to stiffened panels as commonly used in 
fuselage, wing, or empennage structures. 

This report describes in some detail the analysis procedures arid 
rationale for the assumptions used therein. A brief description of the 
sizing methodology is given. Further discussions of the optimization 
procedure and routines are available in References 1-3. Detailed inv-t ruc- 
tions for the use of the code and interpretation of the output from the 
program are given along with several examples in a Separate User's Manual 
(Reference 4). 
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SYMBOLS 


Symbols generally are defined in the text as they fir3t occur. A list of 
the more important ones is given below. 


■ij 


- tangent in-plane stiffnesses of postbuckled plate 


- stiffener spacing 


sk 


rotational stiffness of skin used in stiffener local 
buckling analysis 


- normal pressure load 


u.v.w 


- longitudinal, transverse and normal displacements of a plate 


o o o 
u ,v ,w 


- x,y, and z direction displacements of stiffener shear center 


x,y,z 


coordinate directions 


y o’ z o 


- stiffener shear center coordinates relative to centroid 


y ST’ Z ST 


- stiffener centroid coordinates 


y L ,y R ’ Z R 


coordinates of lines where skin attaches to stiffener in 
torsional/flexural analysis 


th 

vertical distance from panel centroid to i stiffener 
element 


ij 


- plate in-plane stiffnesses 


xi 


th 

- longitudinal extensional stiffness of i n stiffener element 
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^sk * ^xsk “ tan 8 ef it and secant longitudinal extensional stiffness of 
postbuckled plate 


D.^ - plate bending stiffnesses 

D xi - longitudinal bending stiffness of i Stiffener element 

£A X - panel longitudinal extensional stiffness 

El - panel longitudinal bending stiffness 

J 

EA, GJ, C . * stiffener extensional, torsional and warping stiffnesses 


El , El , El - stiffener bending stiffnesses 
yy yz zz 


EAg, EA t - secant and tangent panel extensional stiffnesses 

EIj, , EI T - secant and tangent panel bending stiffnesses 

El - stiffener polar bending stiffness 


W 


j.. constraint 

tn 


K * • % 


- rotational stiffeness of stiffener* equations 21 and 23 


- net rotational stiffness of stiffener, equation 24 


H no no 

yy y 


- transverse moment and normal shear constraints on k 
stiffener 


Mg, M n , Mp - moments due to bow eccentricity, axial load eccentricity and 


pressure 


MSj - margin of safety for failure mode j 
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N x N y N xy 


N xi 


N (k) (k) 

yy V N xy 


N 1 ,M 2 ,N 12 


xi 





R , R 

y 


z 


V 


W 


i 


X 


i 


YCX^ 


£ 


- applied longitudinal, transverse and shear loads 

- longitudinal stress resultant in i^* 1 stiffener element 

- in-plane transverse and shear constraint forces on k 
stiffener 

- average longitudinal, transverse and shear stress resultants 
on po3t-buckled plate 

- equivalent thermal loads 

jL 

- equivalent longitudinal thermal loads in 1 stiffener 
element 

- equivalent transverse thermal loads in i> stiffener element 
and skin 

- standard and modified Euler loads 

- stiffener compression load 

- n" degree Legendre polynomial 

- reaction forces on stiffener at center of rotation 

- strain energy per unit area of postbuckled plate 

th 

- width of i stiffener element, i s 1,5 

fch 

- i design variable 

- objective function 

- laminate strains 
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W a z ~ longitudinal, transverse and normal 3tres3es at 
skin/stiffener interface 

T xy ’ T xz ,T yz " i n "Pl ane and transverse shear stresses at skin/stiffener 
interface 
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STRUCTURAL ANALYSIS 


A flow chart of the principal elements of the analyses is shown in 
Figure 1. Initially skin/stiffener section properties are calculated 
assuming the skin to be unbuckled to obtain a first estimate of the 
longitudinal strain and curvature of the panel. Since the rotational 
restraint provided by the stiffener to the skin i3 a function of the 
3tiffener load, it may be determined only after the panel strain and 
curvature are known. If the applied loads cause buckling of the skin 
(assuming skin buckling is permitted), the postbuckled plate analysis is 
performed and an iterative procedure is needed to determine the correct 
strain and curvature of the panel. The postbuckled plate analysis provides 
the tangent and secant stiffnesses of the buckled skin needed in determin- 
ing the panel strain and curvature. The rotational restraint of the skin 
by the stiffener is reevaluated during the iterative procedure to account 
for redistribution of load between the skin and stiffener. Skin strains at 
various locations are also compared with material allowables and the mem- 
brane strains are checked against prescribed strain limitations. Strength 
checks and local buckling checks are performed for the stiffener element'! 
including the rolling mode in the case of flanged stiffeners. The tor- 
sional/flexural buckling analysis of the stiffeners takes into account the 
effect of the attached skin. A skin/stiffener interface stress analysis is 
performed as a final check at the option of the user. Design constraints 
relative to the skin layup may also be imposed. 

The principal elements of the analyses are discussed in some detail in 
the following sections. 

SECTION PROPERTIES AND PANEL STIFFNESSE S 

The geometry of the stiffener is shown in Figure 2. The width of any 
or all flanges may be set to zero to produce other than I-shaped cross- 
sections. The widths, thicknesses, and ply lay-ups of all five stiffener 
elements may be different. All elements are assumed to be balanced and 
symmetric. It is further assumed that all plate elements have a sufficient 
number of plies so that bending-twisting coupling in the plates may be 
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Figure 1 , Flow CFtort of Analysis Routines 
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neglected, and the plates may be treated a3 specially orthotropic. These 
requirements are normally satisfied by practical composite structures. The 
simplifications allowed by these assumptions greatly increase the computa- 
tional efficiency of the analytical routines. 


Z 



Figure 2. Stiffener Geometry 


The stiffener is treated as an assemblage of plate elements in deter- 
mining the extensional, bending, torsional and warping stiffnesses. Linear 
laminate theory is used to determine the plate inplane stiffnesses, ... 

and bending stiffnesses, D , as well as the equivalent thermal loads, 

T 1 J th 

N . The stress resultants acting in the plane of the i stiffener element 

are given by 


l N ti - Hli - Ml 


( 1 ) 
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I# M assutoed that all of the # loading Is taken by the skin. The 
IbWgi^utfinal stress resultant in each of the stiffener elements is written 
itf form 


N . = ft . 

XI XI 


€ . + 
xs 


T 

ir . , N . 

xip 17 1 


- N 


,T 


XI 


(2) 


W80 tfre longitudinal stress resultant in the skirt becomes 

N xsk = \sk ^xsk *xysk % + N ysk } N xsk 
where the longitudinal extensional stiffness, is 

A x = A 11 " A 12 /A 22 

Similarly for the longitudinal bending stiffness 

\ = D 11 “ °12 /D 22 

The longitudinal extensional and bending stiffnesses of the panel are 
required irt the beam-colurtin analysis which determines the longitudinal 
strains and curvatures in the various elements of the panel (see Figure 2). 
These stiffnesses are: 


FA. = 


El * 


5 

2_i w. A . + b A , 
1 xi s xsk 


(4) 


4 

I w. 


i=l 


A 


(A xi Z i * D xi>J + *5 A x 5 < z 5 + W 5 /12) 


- 2 


+ b s (A xsk Z r;k + »W 


(5) 
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where the z^ distances are measured from the panel centroid to the mid- 
surface of the plate elements. If the skin is buckled, secant or tangent 
values of X ^ must be used in calculating panel stiffnesses. This is dis- 
cussed in more detail in the Postbuckled Plate Analysis Section. 

In the torsional/flexural buckling analysis, extensional (EA), bending 

(El , El , El ), torsional (GJ), and warping (C .) stiffnesses of the 
yy yz zz i 

stiffener are required. In calculating these stiffnesses, the portion of 
the skin which is integral with the attached flange of the stiffener is 
assumed to be part of the stiffener. Elementary beam theory is used to 
determine these stiffnesses as well as the location of the shear center 
about which the warping stiffness is computed. 

BEAM-COLUMN ANALYSIS 

Since the stiffnesses of the buckled skin are functions of the skin 
and stiffener strains and the skin strains depend on the panel stiffnesses, 
an iterative beam-column analysis is required to determine the strains and 
curvatures in the panel. If the skin is not buckled or is raquired to be 
unbuckled, only one pass through the iteration routine occurs in the 
analysis. Otherwise the iteration proceeds until the newly calculated 
longitudinal strain in the skin approximately equals the previously cal- 
culated value. A simple flow chart of the basic operations is shown in 
Figure 3. 

After selecting any initial guess for the skin strain and calculating 
the stiffener properties, the unbuckled panel stiffnesses are computed. 
The Euler load, adjusted for stiffener shear flexibility, is defined as 


P E P E^ 1 ' ^E^^b W 5^ 


( 6 ) 
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Figure S, Flow Chart of Beam -Column Strain Iteration Routine 
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where Agg and w^ refer to the shear stiffness and height of the web, 
respectively, n Is the shape factor, EI^ Is the tangent bending stiffness 
of the skin/stif fener combination, and Pg Is the standard Euler load. 
Equation (6) is obtained from Reference 5, Section 2.17. For stiffeners 
having a free flange (w^ w^), the shape factor, n, is conservatively 

assumed to be 1.0. For blade stiffeners, n is taken as 1.1, a value midway 
between that appropriate for a rectangular cross-section (1.2) and a 
flanged stiffener (1.0). 

The maximum moment at the panel midpoint is the sum of the moments due 
to initial bow, load eccentricity, and pressure 

M = Mg + Mj^ + Mp (7) 


If the magnification parameters 7 and \X are defined as 


y = 


N b /P_ 
X 3 E 


(8) 


H =7Tn/>V2 ( 9 ) 


the moments can be written as 

Mg = -N x b s e /(l-r) (10) 

Mjg = -N x b s A/cos (fi) (11) 

Mp = -pL 2 b s [sec(/i) - l]/4/| 2 (12) 
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(13) 


Mp = ^pL 2 b s Mp)feic(#i) - U/4p 2 (Fixed Ends) 


where 

e s Initial Bow Magnitude at Midspan 

(Positive Inward , Both + And - Are Checked) 

h = Eccentricity of N From Centroid 
(Positive Towards Skin) 

p = Applied Pressure 

(Positive If Internal) 

The curvature at midspan is 


K 

x 




( 1 ) 


where ’El is the secant bending stiffness of the panel. Although the 
tangent and secant stiffnesses vary along the length of the panel due to 
the variation of the moment along the length, the panel is conservatively 
assumed to have constant properties along its length based on the maximum 
mid span moment. If a nonzero value of "e" is specified, separate analyses 
are performed for +e and -e . 


The longitudinal strain in any element of the panel is 


€ . 


xi 


- V®S + Z i *K 


(15) 


in which P is the total equivalent longitudinal load on the panel 
X 


P = b [N -m. , (N + N T , ) + N T , ] + nL, 

x s 1 x xysk y ysk xsk 


( 16 ) 
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T 

where N ^ is the equivalent thermal load in the stiffener. Having an 
improved estimate of the longitudinal strain in the skin from equation 
(15), as well as strains in the stiffener elements, new skin secant and 
tangent stiffnesses can be computed. These stiffnesses are not only 
functions of the skin strain but also of the stiffener element strains as 
reflected through the stiffener rotational restraint of the skin. This 
rotational restraint is discussed in the following section* The iteration 
on skin strain continues until two consecutive values differ by less than 
0.000020 inches/inch. 

ROTATIONAL RESTRAINT OF SKIN 


The skin plate is restrained along it3 long edges by the stiffener. 
The panel ends are assumed to be simply supported. For long wavelengths, 
the rotational restraint is primarily a function of the torsional stiffness 
of the stiffener. For short wavelengths, the web transverse bending 
stiffness dominates the rotational restraint. Intermediate wavelengths 
require interaction of the torsional stiffness of the stiffener and the 
bending stiffness of the web. 

The torsional stiffness of the stiffener may be determined from the 

differential equation of a column subject to a distributed twisting moment, 

M . This equation, obtained from lection 5.5 of Reference 5, is 
a 


c, <j> , 

1 XX XX 


(<17 - El P/EA) 4>, - P (y w, - z V, ) 

o' 'xx 7 o 'xx o 'xx' 


-R (z -z ) + R (y -y ) = M 
y o a z *0 *a a 


(17) 


in which EI q is the polar bending stiffness; p is the compression load 

acting on the stiffener; y , z are the shear center coordinates and y , z 
• o o a a 

are the coordinates of the center of rotation relative to the centroid; and 
Ry and R z are the horizontal and vertical reaction forces at the rotation 
center, respectively (see Figure 2). 
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If the skin is assumed to deform a nt 1 3 y m met r i e a 1 1 y about the rotation 
center and this point lies in tbs sfeln midplane,, tbs vertical reaction R z = 
0 an# £ s-i . (see Figure 2>. The reaction forces are given by, 

""&■ St 


R = El v,° + El 

y zz , mom yz 

o 

w 'xxxx 

* * ; 

it 0, ). 

o v 'xx 

(18) 

R = El v»° + El 

z yz mom yy 

Q 

W f 

f xxxx 

* ^ ■ * 

: 0, ) 
o ^ xx 

(19) 

Assuming that the cross-seat ion F< 

states 

without local 

deformation 



rotation center, as shown in figure H(a), 



Figure 4. Stiffener Rotational Deformation 


Equations ( 17) through (20) can be combined to yield two coupled 

differential aquations in w and 4>. If M , w* and 0 are assumed to vary 

•P 

sinusoidally, the two equations may be solved to yield the torsional 
stiffness 
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= M a /<t> 


2 2 

ff, Iff [Ci + EI zz (z o - Za ) 2 1 + GJ - P (EI^BA - z 2 + z 2 ) 
A 2 I A 2 


-2 

-[EI (2 -z ) - 0 

yz o a \2 


pyj /(ei 


7T 


yy a 2 



( 21 ) 


where A Is the half-wavelength. As A becomes small increases rapidly. 
In the limit the rotation <p vanishes and the stiffener deforms as shown in 

Figure 4(b). The rotational stiffnesses of the stiffener elements can be 
obtained by solving the differential equation 


D. , w, + 

11 'xxxx 


2< D i2 + 


2D 66 I 


u + w. — N w f - 0 (22) 

w 'xxyy 22 'yyyy x 'xx 


with the appropriate boundary conditions (free or elastically restrained) 
at one edge and an applied sinusoidal moment at the opposite edge. The 
stiffness associated with this deformation is denoted 


K 9 = w e /0 


(23) 


and is a function of wavelength, element loads and element stiffnesses. If 
both and Kq are finite, the total deformation is as shown in Figure 
4(c). The net stiffness is approximately that of two springs in series 

K R * 1/U/K0+ 1/K e ) m) 

If the stiffener has no free flanges, cross-sectional deformation is 
negligible and is set equal to . 
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PQSTgMKLED^ PLATE AMELYSIS 


The analysis' predicts the behavior of anisotropic plates loaded- in the 
postbuOkling range if a combination of in-plane biaxial compression, or 
tension*, and shear . Derivations are based: on the 3hear field theory 
developed by W. f. Koiter (Reference 6) for isotropic plates. This theory 
w#s attended to include the case of symmetrically laminated composite 


The tefllrif displacement pattern 
expressed by 


used in the analysis is 


w(x,y) 


sin v (x-my)’ 


inf which A is the half -wavelength of the buckle in the longitudinal ( x)' 
direction and ffl defines the inclination of the noda 1 lines in the presence 
of shear . To extend the analysis into the advanced postbuck ling regime* 
the function W(y) is taken as constant (W s f ) in a center strip of width 
(1 - Of 'b aM remains doubly curved in the edge zones. The bo-ek 1 1 og 
displacement pattern is shown in Figure 5. The Rayl eigh-Ritz method is 
used to determine the unknown parameters. A, m, f and a . For plates of 
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finite aspect ratio, A must be determined so as to produce an integer number 
of half waves along the plate length. The correct number of half waves is 
that which makes the strain energy a minimum. 

Average stress resultants 


The total strain energy per unit surface area, including the torsional 
strain energy of the stiffener, may be written as 




N 2 2U_. N,N- N? N 2 


f - + r + af + I p2 [ < 3 - 2 “> aA u ° 2 + a 22^ 3 j 


A n A n 


A 22 66 


+ ~ n*p[d 2 c {2-a) + | D c 2 /a + U 22 /ot 2 } | +K r F/(a 2 b s ) (26) 


where 


p p p 

P = Ff /(4b ) 


D - b s 2 /X 2 


N* = 47T 2 D 11 /b 2 


= D ij /D ll 


M 3 M 12 + 2/i 66 


c * 1 4/i.^gfn + 2^m - 4/i 26 m 3 + M 2 


C k = a k C/e)rn k 


k = 1,2,3 
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The rotational restraint (R^J is a fun'eti'oh of X and of the Stiffener 
load (see Rotational Restraint of Skin ) . The average stress resultants, 
Hj, Ng, and acting on the postfeuckled panel May fee written in terms of 
the panel edge strain, (y, jig* and and the wave shape parameters D, F, 

m, and ot. 

N 1 = A H € 1 + A 12 f 2 * {A 11 4 \ 2 U -«/2) OF 4 A l2 F/(2a) 

} 

n 2 = A 12 € 1 + A 22 € 2 + (A 12 * A 22 w2) (1 " tt/2) DF + A 2 2 F/(20) (27) 

N 12 = A g g [r i2 - (2-af)fflDF| 

Differentiating the strain energy equation (26), with respect fcfe F ( trii 
and a, for a given value of D* yields three equations from which these 
three parameters can be determined . 

N-jDU-a) + N 2 tA(2 -a) + l/m - m i2 mD(2-a) 4 2K R /(tt 2 b g ) + 

|f [ (3- 2a) aA n D 2 4 A 22 /a 3 ] + i m (D 2 c(2-a) + ±nc 2 A+lJ- 22 /& 3 } = 0 (28) 


N 2 in(2-a) - N 12 (2-a) + | N* {OC^-O) 4 i C^/Ot) = 0 (29) 

N-jD + N 2 (m 2 D + 1/a 2 ) - 2N 12 rrD - |f [ (3-4a) A 1} D 2 - 3A 22 /a 4 ] 

4 ^N*{D 2 C + •|DC 2 /a 2 4 3M 22 /a 4 ] + 4K R /(a 3 b g ) * 0 (30) 


Since equation (30) was obtained fey differentiation with respect tod. 
it Is not valid in the initial postfeuckling region where Ot Maintains a 
value of 1.0 over a range of loads greater than the initial buckling load. 
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Initial buckling of sjtin 

At the onset of buckling ct = 1 and F - 0 . Equations (28) and (29) may 
then be written 


NjD -t- N 2 (m 2 D + 1) -2N 12 mD + (D 2 C + C 2 + + ~ = 0 ( 31 ) 


N 2 m " N 12 + + i C 3 ) ~ 0 


(32) 


For specified values of and N 12 , the quantity m can be determined 
from equation (32) by iteration, for a given value of D (wavelength). The 
above parameters together with a preliminary value of may be substituted 
in equation (31) to yield the corresponding load Since is a 

function of stiffener load a 3 well as wavelength, the correct value of 
again must be found by iteration. 

T,,-' average stress resultant N ^ is calculated for a series of wave - 
lengths. The minimum value of corresponds to the intial buckling load 

N icrr 

Strain s in buckled skiji 

As one of the failure modes considered in the analysis, strains, or 
stresses, in the skin are compared with material allowables. The maximum 
longitudinal membrane strain ( ) in the skin occurs along the nodal lines 

parallel to the stiffener. This strain increases with increasing load. 
The longitudinal membrane strain in the center portion of the skin between 
two stiffeners on the other hand, remains essentially constant and equal to 
its value at initial buckling. Large bending strains exist in this region 
of the skin and, as a result, the total compressive strain in- the center 
portion of the skin may exceed the strain in the skin at the stiffener 



location. In the center region of the skin, the middle surface strains are 
given by 


6 x = « + FD 


€ y = € 2 + FD [in (1- a/2) - a/2] + F/(2a) 


(33) 


y xy = r !2 ~ 2FDm(1 ~ Q/2 


The curvatures in the center portion of the skin vary sinusoidal 1 y in 
the longitudinal direction 


w 

r xx 

- 277DF 

w 


2 

m w. 


yy 

f xx 

w 

» . ■" 

-mw. 


xy 

r 'XX 


(3*0 


The strains in each ply are computed at 'the four potentially critical 

points defined by (x,y)=(0,0) , (A/2. 0), (0, b /2) , and (A/2, b /2) , 

s s 

skin strength margin of safety corresponds to the ply and the location 
where the lowest margin of safety exists according to the maximum strain 
criterion or the Tsai-Hill criterion. 


Buckled skin stiffnesses 


For stability calculations, the stiffness matrix of the buckled skin 
with respect to incremental deformations must be known. In terms of the 
average stress resultants in the skin, and Ng = and the 

strains in the skin at the stiffener locations , ? and €g = T 1 ^, the 

incremental load-strain relationships may be written as: 

a ll a 12 a 16 
a 12 a 22 a 26 

a 16 3 26 a 66 




( 35 ) 



30 that the coefficients of the tangent stiffness matrix are given by 


dN. 

a ij 'j- 1 


i,j = 1/2/6 


( 36 ) 


To determine the coefficients a^, equations (27) are differentiated 


with respect to €j. This yields equations in the form 


where 


i ajl ) 

r 

H 

if 

( a j6 > 

L 


b ll b 12 b 13' 
b 21 b 22 b 23 
b 31 b 32 b 33 


P. 

3 




a. 



) 1 

F ) 

( - d ' 

i 

a €j 

) I 

° ) 


B, 


B, 


B- 


(37) 


j=1,2,6 


Differentiating equations (28) through (30) with respect to gives 


C C C 
11 12 L 13 


c c c 
21 22 ^23 


c c c 
31 u 32 33 


| 

>) 

< 

| aj2 [ 

1 

( a j6 ) 


' d ll d 12 d 13 


d 21 d 22 d 23 


d 31 d 32 d 33 j 


F. 

3 


m. 


a- 


= 0 


(38) 


j = 1,2,6 


Simultaneous solution of equations (37) and (38) yields the tangent stiff- 
ness coefficients a^. 
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Longitudinal tangent and secant stiffnesses of the buckled skin are 
required to calculate section properties for use in the beam column 
analysis. The longitudinal tangent stiffness 

. dN, 

y tan I 

xsk d€^ 

Is determined as follows. Inversion of equation ( 35) yields 

fdC} = [a*] dN 

Since the average stress resultants N ^ and N ^ are assumed constant in the 
post-buckling regime, dNg * dM^ = 0 and hence 

_tan dN-^ ^ 

fi Ksk " m) 


The longitudinal secant stiffness is given by 

^ ° ec _ N 1 “xsk ^xysk ( 'V +N ys k ) 

XS ' k € + fl? - V N 1 ’ 1/A 

1 *- xsk *xysk ysk-*' A 11 

DESIGN STRAIN LIMITATIONS 


(40) 


Structures currently designed with composite materials are limited to 
design ultimate strain levels substantially below the failure strains of 
the basic composite material. The two principal factors establishing these 
limits are (1) stress concentrations associated with cutouts, joints and 
splices for structures loaded in tension, and (2) tolerance for impact 
damage for compression critical structures. Other considerations are! 
transverse cracking in plies perpendicular to the major loading direction , 
temperature and humidity, compatibility with adjacent metallic structures, 
manufacturing defects, and repair ability. 

In POSTOP, tension and compression strain limitations in the fiber 
direction and transverse direction may be specified by the user. These 
limitations are applied at the ply level to the skin membrane strains only. 
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Local bending strains due to pressure loading and postbuckling are not 
included in checks against these limitations. Likewise thermal strains are 
excluded. Membrane strains in the skin resulting from panel bending due to 
eccentricities, pressure, and skin buckling, however, are included and 
added to those caused by in-plane loading. 

STIFFENER STRENGTH 

The strength of the stiffener is based on the maximum strain criterion 
or the Tsai-Hill criterion as applied to the critical element of the stiff- 
ener. This element is the free flange, if one is present, or the free edge 
of a blade stiffener. Longitudinal strain, is computed from equation 
(15). Since the transverse applied load and in-plane shear are zero, the 
transverse and shear strains are 


yi 


« (-A 


12 


f xi - *yi>/ A 22 



T 

where N ^ is the transverse thermal load in the element. If the stiffener 
is a blade, the curvatures of the web about its midplane are set equal to 
zero. If a free flange is present, the longitudinal curvature about its 
midplane is set equal to the panel curvature as given by equation (14). 
All other curvatures are set equal to zero. Strains are calculated at each 
ply and the critical margin of safety is computed according to the speci- 
fied criterion. 

Local bending and twisting strains in the stiffener elements due to 
postbuckling deformations in the skin are currently neglected in the 
analysis . 


STIFFENER LOCAL BUCKLING 


The skin is allowed to buckle locally in POSTOP as discussed pre- 
viously. The wavelength at initial buckling is 2 X it) . As the skin is 

loaded further into the postbuckling range, the critical wavelength 
decreases. POSTOP requires the stiffener to offer positive rotational 
restraint to the skin in all wavelengths shorter than 2A^. 
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The stiffener will therefore always restrain the skin and will not 
tend to drive the skin or to force it to have a stiffener dominated 
wavelength. For the longer wavelengths (A>A it) )» it is assuned that the 
Stiffener may be restrained against local buckling by the skin. For 
panels with normal geometry, this restraint will be advantageous only in 
the case of panels with blade stiffeners. The critical Wavelength for 
local buckling of flanged stiffeners is normally smaller than the critical 
skin buckling wavelength (A<A.^ b ) . Torsional/ flexural buckling will govern 
at the longer wavelengths. The local buckling analysis is performed for 

all half-wavelengths from A* L to A M = L /N - , where N__„ is a user 
defined maximum number of wavelengths. For A^ > A ib the skin rotational 
stiffness, K g . t is conservatively estimated as HDgg/b-, Here Pgg is the 
skin transverse bending stiffness neglecting the additional stiffness due 
to the attached flange. In panels with buckl ed skin, the out-of-pl an e 
deformations of the skin increase the transverse bending stiffness oyer 
that of a flat, unbuckled plate for A^ > A it) . The converatism in the 
estimated rotational stiffnesd, 4D_ 0 /b , may therefore be extreme for 
highly po st buckl ed panels . A summary of the stiffener local buckling 

calculations is shown in the flow chart in Figure 6. 



Figure 6. Flow Chart of Stiffener Local Buckling Analysis 
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The local buckling analysis ,of the stiffener for a given wavelength 
assumes that each stiffener element is uniformly loaded in pure compres- 
sion. This assumption is approximated in the stiffener web of a panel 

loaded in bending due to eccentricities or pressure by taking the load at 
the web midheight as the magnitude of the uniform load. Figure 4(b) shows 
the deformation associated with stiffener local buckling. 

The rotational stiffnesses of the stiffener elements are obtained from 
the solution to equation (22). First the free flange rotational stiff- 
nesses are obtained from the solution for a plate with one edge free and a 
sinusoidal moment applied at the opposite edge. The stiffener web/skin 
junction rotational stiffness, K q as defined in equation (23), is then 
obtained from the solution for a plate restrained by the flange rotational 
stiffness at the far edge and subject to a sinusoidal moment at the near 
edge. The factor by which the actual element loads must be multiplied to 
give neutral stability is found by iteration for each wavelength. Neutral 
stability corresponds to K 0 + K sk = °‘ The margin of safety is equal to 
the lowest load factor minus one, 

ROLLING OF STIFFENERS 


While stiffener stability at short wavelengths is governed by local 
buckling, torsional modes become critical at long wavelengths. If the skin 
is very flexible in bending, as is normally the case in buckled panels, the 
critical long wave modes will be torsional/ flexural with essentially rigid 
cross-sectional rotation. This mode is discussed in the following section. 
If the skin is very rigid in bending, a torsional mode, as shown in Figure 
7 for a rigid skin, could be critical at long wavelengths. This will be 
referred to as the rolling mode. In lieu of a more accurate, coupled 

torsional analysis combining both the rolling and torsional/ flexural modes, 
two simpler analyses are used in POSTOP. 
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Using an assumed displacement field for the stiffener de'fof-Wdtidft itt 
the roiling mode* the total potential energy mSy be expressed arid miftiffliSSSl 
to determine the critical load factoK Here the displacement field asSumed 
id (Reference 7) 


u = -Cz 2 y ^ cos ^ x 
v = Cz 2 Sin x 
w = -2Czy sini ^ x 


(4lj 


in which u; v* and w are the displacements iti the x, and z coordinate 
directions Shown in Figure 7. 

Neglecting behdirtg energy terms in the flange bompahed to the rriembf-arte’ 
terms, the load factor, R, for Stiffener tolling in a given wavelength cart 
Be Writteri as 


R - {t(fi 2 D n 4/5 - 4„ 2 D 12 /3 * 4B 22 /(|) 2 + 16w 2 5 D 66 /3] web 


* (|) 2 w 5 t%l 4 + A x 2 “l 1 ?/ »5 tfl x9 V 5 + M kl w i + W 2> 


( 42 ) 


The lowest value of the load factor id found fori wivilerigthd ranging 
frSrii the panel length to the stiffeneh Spacing. fhl margin 8f Safety is 
equal to the lowest load factor minus one . 
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TORSIONAL/FLEXURAL BUCKLING OF STIFFENERS 


In the analysis it is assumed that the stiffeners are cocured or 
bonded to the skin and the portion of the sKin that is integral with the 
stiffener flange is considered part of the 3tiffener» The remaining skin 
on each side of the stiffener is assumed to be detached and replaced by a 
set of equivalent stress resultants. These stress resultants (and couples) 
are considered to act as external loads on the stiffener, as shown in 
Figure 8. The buckling mode is assumed to be periodic along the width 
direction with one period ’p' including one or more stringers (see Figure 
9). 




Figure 9, Transverse Mode Shape 
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la determining the constraint forces , membrane and flexural effects 
are assumed to be uncoupled and treated independently. A section of skin 


th 


. th 


between the k and (k+1) stiffener of width b^ a b g - b g ^. is analyzed 


Calculation of membrane forces acting on stiffener 


— (k) — Ck) *- 
If U and v represent the displacements of the skin in the x and 

y directions measured with respect to the prebuckling configuration, 
inplane equilibrium yields the following two governing differential 
equations. 


a 77 W , a 77 <k) . » - (k) 

A 11 u, xx A 6€ u 'yy + (A 12 A 66 } v 'xy 


* A 12 + h 66 ) 
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(k) 

'xy 


+ A 


66 
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'XX 


+ A 


22 'yy 


fc 0 


By letting 


u = {jW (y) cos nx 

~ ( k) _ V (k) J jyj 


(44) 


in which n = n7T/L, and substituting equation (44) into equatioh (43). 

(k) 

the following uncoupled fourth order differential equation in U i& 
obtained . 


u 1 


-2 
- n 


|A ll /A 66- 2A 12 /A 22> U 


.!<«, -4 
+ n 


( a h/a 22 )u 


(k) 


a 0 


(45) 


The solution to equation (45) can be put in the form 


4 

U (k) = 2 A^ k) f.(y) (46) 

j=l 3 3 
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where f . are known functions of y and A 
J J 


(k) 


are constants of integration. 


Substituting equation (46) into the first of equations (43) and integrat- 

(k) 

ing yields an expression for V which may be written in the form 


V 


(k) 


= I 


j=l 


A j 


(k) 


9j (y) 


(47) 


The integration constants may be expressed in terms of the buckling 

L U tU 

displacements of the k and (k + 1) stiffeners by imposing conditions of 
continuity between skin and stringers at the left (L) and right (R) 
attachment lines shown in Figure 8. Denoting the longitudinal displacement 
of the stiffener during buckling by u°, the displacements of the shear 
center in the y and z directions by v° and w° , respectively, and the twist 
per unit length by $ , the displacements of the stringer at an arbitrary 
point q of the cross-section may be expressed as 



. p 

V i 


_ . p . 

Z q 


(W. 


- W ) 

q X 


V q = V ° ‘ 
w q " w ° + 


(48) 


Ihe coordinates of the shear center (y Q , z q ) and those of point q (y q , 
Zq) are measured with respect to the centroid. The last term in the first 
of equations (48) represent the displacement due to warping of the cross- 
section at q. To satisfy inplane continuity between skin and stiffeners 
at the attachment lines, the following conditions are prescribed. 



( 49 ) 
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By substituting equations ( 46) and (47) into Equations (44) and using 
the results together with equation (48) in the boundary conditions (49), 
one may solve for the integration constants in the form 



(50) 


where the terms in the 4x4 matrices [R] and [L3 are functions of f j(?) 

and gj(y) at y : 0 and y s b^» respectively. 

If there are N stiffeners within one period, there will be N equations 
til 

of this type. The N equations, however, becomes 


(N) 


u 


u 


,(i) 


CCS ;iX : l ft) 


V, 

w, 
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'x 

o 


(51) 


(k) 


(k) 


th 


The constraint forces N and N acting on the k stringer may nOw 

yy xy 

be determined in terms of the displacements and twist of that Stringer and 
the adjacent ones, or 


(10 o<K> , <*> 

= s. u ° + v° + s„ +s_ a 


yyR R^'x ' 'xx 


R^ 'xx R 4 'Xk 


n (k+l) (k+1) (k+1) . (k+1) 

+ s l. u 'x + V'« + V'» + V'** 

1 2 3 4 

/uy ^ n (k-D ^ n (k-l) _ (k-1) _ (k-l) 

N' k l *= S„ u,° + S_, V,° 4- S n w,° + S_, 0/ 

yyL R-^ x R^ 'xx xx R^ 'xx 


_ 0 <K> - 0 (K) _ 0 (k) .. (k) 

+ S L 1 u, x + s t 2 v 'xx + s l 3 m 'xx + s l/'*x 
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Ik) o (k) o (k) o (k) (k) 

n x?r * v + + V'" + V* 


(k+1) 


(k+1) 


(k+1) 


+ T r u° + T t v ° + T t w ° + T, 0, 


(k+1) 


(52) 


l L 2 "x ’ "L 3 ” f x ‘ *L 4 ^x 


(k) - o (k_1) - 0 (k-D _ rt (k-l) _ (k-1) 

N = T u + T v ° +Tw^ + T d) 

XVL R. P V 'v + A t> + J-n 

I X 

(k) > 


1*2 'x ' ^"'x ' ‘R 4 "x 


o (k) = o (k) - o (k) 

- * m .. v i m . „ 


+ + T L 2 V 'x + T L 3 “'x + T l/’> 


where the terms S_ , jL , T_ , and T_ are functions of the matrix 

R j R j R J . J _ 

coefficients R, , . The terms S. , S. , T. , T. are functions of the matrix 

J j j j j 

coefficients Ljj . 

Effect of skin bending on stiffener buckling 


The governing differential equation for plate bending in the presence 
of membrane forces (R , N , N ) may be written 

a A j Jr A J 


D ll w 'xxxx + ®16"'=c^ + 2 < D 12 +2D 66> + 4D 26 u 'x^ 


+ D-_w, ^ - N w,< k > - 2N w/ k ^- N w, ^ - p^ k * = 0 
22 YYYY XX xx xy xy yy yy * 


where w is the normal displacement of the skin panel with respect to the 

pre-buckling configuration and p is the normal pressure. 

In determining the flexural effects of the attached skin, the stress 

resultants N and N are treated as constant loads. In order to sim- 
xx yy 

plify the analysis, membrane shear R - is not included and the bending- 

xy 

twisting stiffness terms and are assumed to be zero. In the 

absence of normal loading, the governing differential equation reduces to 


D u w, 


< k > + 2D,w, < k > + D„,w. < kl -S w, (k >-N w, (k) = 0 


xxxx 


J *v r 1 / i-v w / tv rv / 

3 xxyy 22 yyyy xx 'xx yy 'yy 


(53) 


where 


°3 = °12 + 


2D, 


66 
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Representing the normal displacement as 


and defining 


_(k) (k) 

w W (y) sin nx 


_2 

g = [D 3 + N /{2n )}/^ 22 

i-P 2 = (D u + V 5 2 )/(I> 22 9 2 ) 


equation (53) may be written in the following form. 
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,iv 


(k) 


- 2g n W 


ii (k) 2 2 

+ g (i -in 


-4 

a 


(k) 

W = 0 


( 54 ) 


The solution of the differential equation (55) depends on the values 

_2 

U and g. The general form of the solution may be written as 


W (k) = S C, (k) F. (y) (56) 

j=l 3 J 


The integration constants C 

th 


(k) 


displacements of the k and (k 


are solved in terms of the buckling 
th 

+ 1) 1 stiffeners by requiring the 


following conditions to be satisfied 
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Using the third of equations (48) in conjunction with equations (54) 

th 

and (56) one obtains the solution for the integration constants of the k 
stiffener in the form 



where the terms in the 4x4 matrices [IT] and [C] are functions of the 
F^(y). For the N fch stiffener 



The moments and shear forces acting on the k stiffener can now be 

th 

obtained in terms of the normal displacements and twists of the k 
stiffener and the adjacent stiffeners. 



The terms B n , B„ , R_ , and R D are functions of the matrix coefficients 
J J "j J 

17... The terms B, , B. , R, , and R, are functions of the matrix 
1J _ L j L j L j L j 

coefficients j . 
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Formulation of stability matrix 


Differential equations are obtained by considering equilibrium of 
forces and moments acting on a small element of the stiffener. Included 
in these forces and moments are those applied along the left (L) and right 
(RJ attachment lines by the skin during buckling. 

The moments, M and M , and the torque M , may be expressed in terms 

y z ^ 

of the stiffener displacements and twist 

M = - El v,° - El w, ° 
y yz xx yy xx 


M = - El v, ° - El w, ° 
z zz xx yz xx 


( 61 ) 


M = GJ 0, - C. <p, 

X X 1^'xxx 

The following four stability equations are derived by considering 
equilibrium of forces and moments 
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The expressions for the constraint forces and moments, equations (52) 
and (60), may be substituted in equations (62) to yield a new set of 
equations for the k th 3tiffener. This new set of equations will involve 
the buckling displacements and twist of the adjacent stiffeners as well as 
those of the stiffener under consideration. 

Buckling displacements are assumed in the form 

_(k) (k) 

u = U* (cos nx 

Q (k) (k) 

v = V* sin nx 

(k) (k) _ (63) 

w = W* sin nx 

4> m sinra 


Equations (62) may now be written in the form 



(64) 


In the above [AA3 is of order 4N x 4N and [BB] is of order 4N x 3N. 
By eliminating the column submatrix { U* } a standard eigenvalue problem 
of order 3N x 3N is obtained. 



(65) 


The lowest eigenvalue, F , corresponds to the stiffener buckling 
load. 
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skin/stiffener interface stresses 

When the panel is subject to internal pressure loads or has a buckled 

skin and when the panel stiffeners are bonded or cocured to the skin* the 

critical failure mode may involve Separation of the skin and stiffener. 

Prior to separation, the interface between the skin and stiffener attains 

a complex stress distribution. This distribution of stresses may include 

longitudinal and transverse direct Stresses* ff and 0 , normal direct 

x y 

stress, , inplane shear stress, T , and transverse shear stresses* T xz 

and t . 
yz 

The interface stress analysis in POSTOP is an attempt to provide the 
designer with a means of minimizing the tendency toward skin/stiffener 
separation. Prevention of this mode from being critical , rather than 
prediction of the separation failure load, is the objective of the current 
analysis. The computational efficiency required of any analysis in a 
practical sizing code is a strong constraint on the types of analyses 

which may be included in Such codes. This requirement excluded from 

consideration any highly non-linear separation failure analysis. 

The current interface stress analysis provides an efficient, 

closed-form solution for the interface stresses based on the assumption of 
elastic linear material behavior and small deflections. Margins of safety 
are obtained with the Tsai-Hill’s criterion as applied to a three 

dimensional stress state. 

The physical model for which the solution is developed is shown in 

Figure 10. The model consists of flange and skin plates and an interface 
layer. The attached flanges are assumed to be of equal width b and equal 
uniform thickness, t^.. The remainder of the stiffener i3 replaced by a 
vertical and a rotational spring in the web midplane having stiffnesses k z 
and k r , respectively. The skin between stiffener flange edges is replaced 
by equivalent moments and shear forces in the buckled skin. These stress 
resultants are distributed sinusoidally over the buckled skin wavelength, 
A. A transverse stiffness k^, representing the transverse membrane 
stiffness of the removed skin, acts on the edges of the remaining skin 
plate. Longitudinal inplane loads are present in the flange and skin and 
a transverse inplane load in the skin. Inplane shear is currently not 
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Included in the model. The effect of in plane shear is, however, included 
in determining the buckled skin equivalent moments, M, and shear forces, V, 
applied to the skin portion of the model as shown in Figure 10. 


INTERFACE 



Figure 10. Skin/S^iffener Interface Model 


Independent solutions have been developed for the antisymmetric 

conditions present in postbuckled biaxially loaded panels and for the 

symmetric condition that results from pressure loading. The antisymmetric 

case closely approximates the deformation of postbuckled panels loaded 

with biaxial plus shear loads if the shear load is of a secondary nature 

as previously assumed and the flange width, b^., is small relative to the 

wavelength. Since the solutions are linear, superposition is used to 

obtain total stresses if required. The interface normal stress, O z% and 

transverse shear stresses, T and T , are obtained directly. The 

x z yz 

remaining interface stresses are obtained from the skin inplane strains 
and the requirement of skin-interface compatibility. 

The solution procedure for the antisymmetric case is summarized below. 
The procedure for the symmetric solution i3 similar. 

If the flange and skin are orthotropic and loaded antisymmetrically 
about the web midplane, they will deform antisymmetrically. The applied 
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moment and shear may be represented in series form as 


M(|) = sin m7T $ 

m 

V($) =%V m sin rfi7T| 
m 


in which m = 1, 2, . . . and | = x/X * The Skirt displacements in the x, y, 
and i directions may be written in Series form as 

u = £ Z % p n+1 W cos m7r ^ 

m n 

V= II V mn p n (7j) sin (67) 

m n 

“ ' Z Z "Wl W» sln 

m n 


in which n = 0, 2, 4, ... and fj = y/b f . The P ft ( 77 ) are Legendre 
polynomials. The interface stresses acting bn the upper surface of the 

skin may be written in similar fori 

T xz ZI p* W» «* “5 
m n 


T yz 


>Z 

m 


p,__ P„(f|) sin 
ymn n ' 


(68) 



%in p n +i(’7) sin mir% 


A weighted residual procedure is used to establish the governing 
equations for the irtplane and transverse displacements of art orthotropic 
plate subject td surface normal and tangential loads Representing 
interface stresses. iiaxiSl inplane loads are included in the governing 
equations. Substitution of the series expressions for the applied edge 
moments and shears, the skin displacements, and the interface stresses 
into the governing equations yields a Set of linear algebraic equations 



for the displacement coefficients as functions of the applied edge load 
and unknown surface load cpeff icients. In matrix form these equations are 



[g] {«} = 



M + M 


(69) 


in which the terms of the coefficient matrices are functions of the plate 
stiffnesses and dimensions. Due to the properties of the Legendre 
polynomials chosen as the shape functions, these terms may be rapidly 
computed. The {fi} and {V} are functions of the applied edge loads. If N 
shape functions are included in the series expressions, equations (69) 
contain 3N equations. These equations may be developed independently for 
each value of m as required by the series expressions of the edge loads, 
equations (66). Similarly, a set of equations may be developed for the 
flange displacements u^., v^., w^.. 



(70) 


The two plates are forced to have compatible deformations through the 
continuity conditions 


u f = u -3s (t f + t s ) dw/dx + r z t a /G - (h t f t a /E)acr z /dx 

V f = V -35 (t f + t s ) dw/dy + T yz t a /G - (3s t f t a /E)d<7 z /dy (71) 


w f = w + (t a /E)(j z 


in which t „ , t , and t are the flange, skin, and interface thicknesses, 
is a 

respectively, and E and G are the elastic moduli of the interface layer. 
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If equations (69) and (70) are solved for the displacement coefficients 
and substituted, along wih equations (67) and (68), into the continuity 
cohditions, equations (71), a set of 3N linear algebraic equations for the 
interface stress coefficients results 



These equations may be developed and solved for each value of m in the 
edge load series yielding the interface stress coefficients. The 
interface stresses are then obtained from equations (68). 

The interface stresses are evaluated at a number of points on a grid 
defined by the program user. Typically points along the flange edge will 
be critical although locations near the web may be critical in some cases. 
Some examples are given in the User's Manual, Reference 4. 

Margins of safety at the grid points are computed according to the 
Tsai-Hill criterion. This criterion may be written an 



in which the bar superscript denotes allowable stresses. Different 
allowables may be specified for tensile and compressive direct stresses. 
The allowable stresses appropriate to the sense of the computed stresses 
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are used in computing the values of S at the grid points. The margin of 
safety in thi3 mode is 1/\/$ -1. 

The in plane interface stresses used in equation (73) are 


E 


Or = 
x 


1 - V 


- < £ x * -V 


01 . = 


E 


(i*„ + €„) 


y ! _ „2 y 


r = Gy 
xy xy 


in which € x is the longitudinal strain at the interface midsurface. The 

strains 6 and y are conservatively set equal to the skin transverse and 
y xy 

shear strains, respectively. 

SKIN LAYUP DESIGN CONSTRAINTS 


It is often desirable to place constraints on the proportions of skin 
material oriented in various directions. For example, the "soft skin” 
concept has been used as a means of enhancing the damage tolerance of 
stiffened panels. In this case the skin is required to have a high 
proportion, say 60 percent or higher, of its total material oriented in the 
+45-degree directions. Without a lower bound on the proportion of 
+45-degree material, the sizing procedure would most likely result in less 
than 60 percent +45-*degree material in the skin. 

To achieve designs such as the soft skin concept, lower bounds may be 
placed on the proportions of skin material oriented in three general 
directions or zones as shown in Figure 11, Zone 1, the longitudinal zone, 
is defined as the x-direction (stiffener direction), plus or minus a small 
angle 9. The angle d may be zero or any small positive number of degrees. 
Zone B, the transverse zone, is defined as the y-direction plus or minus 9. 
Zone 2, the intermediate zone, includes directions excluded from Zones 1 
and 3. To achieve the soft skin design described above, lower bounds of 
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0.0, 0.6 and 0.0 would be set for Zones 1, 2, and 3, respectively, and 

d would be set to zero degrees. The purpose of the angle Q is to allow 
slightly off-axis, say +5-degree, material to be used as the basic 
longitudinal material rather than only zero-degree material. Transverse 
material would, in this case , be +(85/95 )-degree material and Q would be 
set at 5-degrees. If no skin layup design constraints are desired, all 
lower bounds are simply set to zero. 



Figure 1 1 , Skin Layup Design Constraint Zones for 
Skin Material Oriental ions 


MULTIPLE LOAD CASES 

P0ST0P has the capability to analyze or size panels with up to five 
separate load cases. In the sizing mode, a minimum-weight panel is 
obtained that satisfies all design requirements for all the load cases. 
The margins of safety for each failure mode in each load case are computed 
and the critical margins of safety are used to formulate design con- 
straints during sizing. 

Each load case has its own inplane and normal loads, eccentricities, 
temperature, required panel stiffnesses, skin buckling requirement, 
material properties and allowables and strain limitations. Stiffener 

dimensions and spacing, laminate configurations, and skin layup design 
constraints are common to all load cases. 
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SIZING 


The determination of values for the panel variables so that all mar- 
gins of safety are within prescribed bounds and the panel has the lowest 
possible weight is referred to as sizing. POSTOP uses the program CONMIN 
(Reference 2) to determine how the variables should be changed from an 
initial design. The program COPES (Reference 1) is used to impose limits 
on the variables or margins to safety, to link variables, and to provide 
design aids such as the sensitivity study option or the approximate opti- 
mization techniques described in this reference. COPES also provides the 
ability to size a panel to meet objectives other than minimum weight. For 
example, a panel with a given weight could be sized to provide a maximum 
positive margin of safety in a particular failure mode. The margins of 
safety are computed in all cases by the analysis routines previously dis- 
cussed. 

OPTIMIZATION PROBLEM STATEMENT 


The standard statement of an optimization problem is: Find the values 
for the design variables X^ such that some objective function Y(X^) is 
minimized subject to the constraints Gj(X^) < 0 and (Lower Bound < X^ < 
(Upper Bound) 

Normally the objective function is the panel weight per unit plan form 
area. Alternatively any margin computed by the analysis routines may be 
chosen as the objective function in which case Y(X^) = -MSj. The design 
variables may include any of the stiffener flange widths, the stiffener 
height, the stiffener spacing and thicknesses of the various laminae 
defining the stiffener and skin laminates. Panel length is normally held 
constant but may be chosen as a variable if desired. 

COPES allows upper and lower bounds to be placed on the variables or 
any function of the variables. Variable constraints are treated as linear 
side constraints. Function constraints may be linear or nonlinear. In 
POSTOP the user defines bounds on all variables and margins. Panel weight 
may also be bounded in a margin maximization sizing. The constraints 
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Gj(Xj) may be expressed in standard form as 


Lower Limit - Value = G, (X.) £ 0 

I ' 

Value - Upper Limit = (X.) s» 0 

COPES internally formulates the constraints in this way. Normally 

each margin of safety will be required to be greater than zero, or some 
minimum value, and will not have an upper limit imposed. In this way 
different margins of safety may be required to have different minimum 
values. Variables must be constrained to have positive values. These 
side constraints are imposed directly by limiting possible values that the 
variables can assume. 

OPTIMIZATION PROCEDURE 


The nonlinear mathematical programming procedure used in CONMIN is a 
modification of the method of feasible directions (Reference 3)- This 
procedure is well documented in Reference 2 and only a brief overview is 
given here. 

A starting design is defined which may be in the feasible or infeas- 
ible region but. should be "reasonable.'’ Gradients to Y(X^) and the 
critical Gj(X^) are computed, POSTOP requires CONMIN to compute the 
gradients by finite differences since most of the margins of safety are 
computed by algorithms rather than explicit functions. Assuming the 
starting design to be feasible, it is modified using the conjugate direc- 
tions method until a constraint is encountered. Further modifications are 
made in a direction that will decrease Y(X^) but not violate the Gj(X^) . 
When no further modifications are possible the optimum design has been 
found. A global optimum cannot be assured. Different starting designs 
may be used to achieve some confidence as to the global or local nature of 
the optimum designs. In most cases starting designs that are well within 
the feasible region allow the optimization procedure to produce designs 
that are very close to the apparent global design. Experience with POSTOP 
has shown that local minimum values of panel weight may be found that are 
approximately equal, but the corresponding designs may be significantly 
different . 
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